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Within the reweighting approach, one has the freedom to choose the Monte Carlo action so that it provides a 
good overlap with the finite-/! measure but remains simple to simulate. We explore several choices of action in 
the regime of small /j.. Simulating with a finite isospin chemical potential — jj, gives a better overlap than the 
standard choice fJ. = 0, with no computational overhead. 
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1. INTRODUCTION 

Because heavy ion collision experiments at 
RHIC correspond to a chemical potential /i as 
small as lOMeV, lattice QCD simulations at small 
fi have attracted renewed interest. To address the 
sign problem caused by the complex fermion de- 
terminant, three approaches have been proposed. 
(i) Analytically continuing results obtained at 
imaginary chemical potential (m) Measuring 
the coefficients of the Taylor expansion of each 
observable about /i = H,^; (Hi) Measuring 
the observables from a reweighted ensemble 
(i) trades the sign problem for that of the an- 
alytic continuation; the approach works well for 
smooth functions of /i Q. We will not discuss 
it further here, (ii) has no sign problem, but 
rapidly becomes complicated as the Taylor order 
increases, {in) is enjoying a successful revival in 
a two-parameter version, where results at fi) 
are obtained by reweighting the (/3o, fJ-o — 0) sim- 
ulation The accessible range of chemical 
potentials and volumes is still limited by the sign 
problem. But also, ensuring a sufficient overlap 
between the desired ensemble at (/?,//) and the 
Monte Carlo ensemble at (/3o,Mo — 0) is notori- 
ously difficult. The acceptable shift of parameters 
normally decreases like 1 / volume. 

It seems useful then to consider the sampling 
of an approximate finite-/Lt action. Reweighting to 



the exact action will be mild and safe. The dif- 
ficulty is to capture most of the exact measure, 
while avoiding the computation of the exact, com- 
plex determinant at each Monte Carlo step. We 
explore three variants of this strategy below. 

2. SIMULATIONS AT SMALL n 

The QCD partition function with chemical po- 
tential, here for staggered quarks, is given by 



Z(m) 



dU det M{^l)— exp(-S'G) 



(1) 



with Sq the gluon action and M{fi) the Dirac 
matrix. Writing det M(/i) = | det M(/i)| exp(z^) 
and defining 6 = ^0, Eq.(|^) becomes 

Z{^) ^ ( dU I det M(/i) I ^e^» exp(-S'G). (2) 



To allow Monte Carlo sampling over a positive 
measure, the expectation value of an observable 
0(/i) is recast into 



( 0(m) ) = 



( e'V/ )z 



(3) 



{ ■■)^ means that the average is taken over con- 
figurations sampled from the partition function Z 



Z 



df/ |detM(Ai)|T exp(-5G) /(C/) (4) 
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where /(C/) is any positive functional of the gauge 
field, e*^// is then the correction (reweighting) 
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factor applied to each configuration. The stan- 
dard strategy corresponds to Z = = 0), i.e. 



/ 



detM(^ = 0) 



det M (fi) 



(5) 



One would like to find an optimal choice for f{U). 
Setting aside the algorithmic issue of sampling 
a given f{U), the optimal choice is that which 
maximizes the statistical accuracy on ( 0(/i) ) for 
an ensemble of N uncorrelated configurations. If 
0(/x) is not specified, it is reasonable to minimize 
the fiuctuations in the reweighting factor, i.e. the 
relative error on the denominator of Eq.(|^), since 
it propagates to all observables. From the central 
limit theorem, that error is, for A'' large 



(icOS0)2)/(icOS^?)2 



(6) 



This expression is minimized for 

/ = |cos0| (7) 

which reduces the denominator of Eq.(|^) to the 



average sign ( sgncosf^ 
be desirable to sample from 



Therefore, it would 



■^opt 



dU \ det M{n)\— I cos 6*1 exp(-S'G) (8) 



The drawback of this sampling choice is that 
the phase 9 of det M(/i) must be evaluated at each 
Monte Carlo step, at a computational cost cx 
for a lattice of size x Nt. Thus, this strat- 
egy seems prohibitively inefficient. We consider 
three alternatives, which can be implemented at 
reduced computing cost and are closer to the op- 
timal choice Eq.(^ than the standard Eq.Q. 
• Method A: / = 1 

The sampling measure is cx det~ [Mt(^)Af(^)]. 
Since det M^(/i) = det M(—/i), this amounts to 
simulating a finite isospin chemical potential Q 
fii = fi. Compared with the standard sampling 

measure oc det~ = 0), this method cap- 

tures, with no computational overhead, the fiuc- 
tuations of the magnitude of the determinant 
which may account for a good part of the rele- 
vant physics at small /i. Indeed, the critical line 
Tc{fJ.) in the isospin-/i case seems to have a similar 
curvature as in the usual isoscalar-/x case M. 



fjU) 
1 



fi^O.l 0.15 0.20 0.25 0.30 
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0.46 


2.21 


15.8 


0.008 


0.05 


0.46 


14.4 


14.2 



Table 1 

Standard deviation (numerator Eq.(^)) of the 
reweighting factor cos 6*// as a function of /i, for 
various sampling choices f(U). This number is 
proportional to the relative error on a smoothly 
varying observable, for a given statistics. The top 
line is the standard choice. 

To verify the improvement of the overlap with 
the desired measure, we generated ensembles of 
200 configurations (/3 = 4.8, = 2, ma = 
0.025, 4"^ lattice), using the standard method = 
0) and using method A, i.e. at isospin /z/ ^ 0. 
Reweighting was then performed on each ensem- 
ble to obtain results at isoscalar chemical poten- 
tial II — iJLj. Fluctuations of the reweighting fac- 
tor were measured by the numerator of Eq.(^). 
Table 1 shows that, for small /i, the standard 
method (Irst line) gives larger fluctuations, i.e. 
a poorer sampling, than method A (2nd line). 

Still, in either method the factor cos 9 must be 
computed for each configuration in the ensemble. 
For small /x, it is possible to approximate 9 by its 
truncated Taylor expansion: 



-Im Tr log Af (^) ^9 + 0{ii^), 



^ u Im Tr M-^M 
4 

dM 



/j=0 



(9) 



where M = — — . The advantage is that the trace 
can then be estimated using n noise vectors rji as 



Tr M-'^M 



1 



T]l M-^M r]i 



(10) 



This approach has been used in Q], with n = 10, 
on a 16^ X 4 lattice. To test the quality of this ap- 
proximation, we measured the exact phase 9 and 
its linearized approximation 9 on each configura- 
tion of our isospin-// gauge ensembles. The scat- 
ter plot Fig. 1 shows that the correlation between 
9 and 9 is excellent at /i = 0.1 (i.e. /i/T = 0.4), 
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Figure 1 . The linearized phase 9 of the determi- 
nant vs the exact phase 0, for 3 values of /i. 

and very poor at fi — 0.3. One may expect this 
correlation to persist for larger volumes. There- 
fore, using the linearized approximation 6 up to 
/i/T ~ 0.4, as advocated in appears justified. 
On the other hand, the stochastic estimator of 9 



Eq.(lO) was extremely noisy and required several 
thousand noise vectors. 

• Method B: / = e 

To approach the optimal / = | cos 6*1 at small /i 
without computing higher derivatives in the Tay- 
lor expansion Eq.(^, one can use cos 9 w 1 — 
~ exp(-i^2). rpj^g additional term Sg = \9'^ 
in the action can be included in an R-type algo- 
rithm |10 , where the Trlog is estimated at each 
step using noise vectors. Here, two uncorrelated 
noise vectors are necessary to estimate . Main- 
taining the St^ accuracy of the R-algorithm in the 
stepsize seems less obvious, but may be feasible. 

To assess the advantages of method B over 
method A, we compare in Table 1 the standard 
deviation of the reweighting factor cos9/f. A 
large reduction is seen at small fi. 

• Method C: / = |cos^| 

To approach |cos0| even better at small /z, one 
may consider the choice / = | cos^|. As above, 9 
can be estimated via Eq. ([l0|) . Then, an unbiased 
estimator of cos 9 can be formed by a stochastic 
Taylor expansion as in which can be used in 
the Monte Carlo update. As long as cos 9 remains 



positive, i.e. for small fi, an important gain is 
possible as seen in Tabic 1. 

3. CONCLUSIONS 

We have considered three choices of Monte 
Carlo actions which can be used to obtain finite- 
results after reweighting. Compared to the stan- 
dard choice of simulating at = 0, they provide 
better statistical accuracy at small {f^/T < 0.5), 
because the reweighting factor fluctuates less. In 
the case of Method A (isospin chemical poten- 
tial fii = fj,), this improvement is achieved with 
no computational overhead. Moreover, results at 
finite fii come for free. 
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